Spin excitations in Ko^Fei.ggSeg superconductor as studied by Mossbauer spectroscopy 
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Mossbauer spectroscopy was used to probe the site specific information of the Ko.84Fei.ggSe 2 super- 
conductor. A spin excitation gap, AE «5.5 meV, is observed by analyzing the temperature dependence 
of the hyperfine magnetic field (HMF) at the iron site within the spin wave theory. Using a simple model 
suggested in the literature, the temperature dependence of the HMF is well reproduced, suggesting that, 
below room temperature, the alkali metal intercalated iron-selenide superconductors can be regarded as 
ferromagnetically coupled spin blocks that interact with each other antiferromagnetically to form the ob- 
served checkerboard-like magnetic structure. 

PACS numbers: 76.80.+y, 74.10,+v 



INTRODUCTION 

There has been a renewed interest in the iron-based 
superconductors since the discovery of superconductivity 
(SC) at about 30 K in the A , 8+s Fe 1 , 6+p Se 2 (A = K, Rb, 



Cs or Tl/K) 1JJ-|4|] compounds due to their unprecedented 
physical properties, such as the coexistence of high temper- 
ature SC with strong antiferromagnetic (AFM) order Q3l 0- 
[H]. However, whether SC and AFM order coexist micro- 
scopically or SC only occurs in the non-magnetic phase 
is still highly debated since some reports support the co- 
existence picture fl IsM Sfl and others favor the phase sep- 
aration scenario Il9l4l2ll. Local probe techniques such as 
muon-spin relaxation/rotation (mSR) H and Mossbauer 
spectroscopy (MS) ifTll [t3l4l5ll have shown that a two- 
component picture is inescapable to describe the system 
correctly, namely, all samples are phase-separated into a 
major AFM phase and a minor paramagnetic (PM) phase. 
Nuclear magnetic resonance (NMR) II16II and fiSR llTTl ex- 
periments reveal that the PM phase becomes superconduct- 
ing below T c . However, the scenario that only the PM 
phase becomes superconducting alone can not explain all 
the above mentioned experiments. Thus further studies on 
these compounds are desired to settle the debate. In this 
case, an investigation on the local magnetic property is 
helpful in understanding the correlation between SC and 
magnetic ordering of these systems. 

MS has been proved to be a very useful tool to probe 
local specific information of the iron-based superconduc- 
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THEH. Especially, when possible coexistence 
of magnetic order and SC presents in the same sample, a 
MS study might reveal rich information. So far, only a 
few work using MS to study these materials [[Til [l3i [LUl 
have been reported. A detailed study focusing on the tem- 
perature dependence of the local magnetic field at the iron 
site near the superconducting transition temperature is still 
missing, which might hold the key to understand the inter- 
esting interplay between AFM order and SC. Therefore, in 
the present work, MS was used to study the magnetic struc- 
ture and temperature dependence of the hyperfine magnetic 



field (HMF) at the iron nucleus of K .84Fe 1 .g g Se 2 single 
crystals. The results provide evidence that a spin excitation 
gap opens up before entering the SC state. Using a sim- 
ple spin model, we show that the ferromagnetically cou- 
pled (FMC) four spins can be viewed as a net spin, which 
couples antiferromagnetically with each other to form the 
checkerboard-like AFM structure. 



EXPERIMENTS 



Single crystals of potassium intercalated iron-selenides 
of nominal composition Ko.sFe 2 Se 2 were grown by the self 
melting method similar to previous reports [QI 0]. Sto- 
ichiometry of high purity K pieces, Fe and Se powders 
were mixed and put in a sealed quartz tube. The samples 
were heated to 1273 K slowly, kept for 2h, cooled down 
to 973 K at the rate of 5 K/h and then furnace cooled to 
room temperature by shutting down the furnace. The re- 
sulting plate-like crystals with a shiny surface are of a size 
up to 6 x 4 x 2 mm 3 . The actual composition is determined 
to be Ko.84Fei.ggSe 2 by energy dispersive X-ray spectrum 
(EDXS). 

Single crystal x-ray diffraction (XRD) measurements 
were performed on a Philips X'pert diffractometer with Cu 
K a radiation. AC susceptibility measurements were car- 
ried out through a commercial (Quantum Design) super- 
conducting quantum interference device (SQUID) magne- 
tometer. Transmission Mossbauer spectra were recorded 
using a conventional constant acceleration spectrometer 
with a 7-ray source of 25 mCi 57 Co in palladium matrix 
moving at room temperature. The absorber was kept static 
in a temperature-controllable cryostat filled with helium 
gas. The velocity of the spectrometer was calibrated with 
a-Fe at room temperature and all the isomer shift quoted 
in this work are relative to that of the a-Fe. 
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RESULTS AND DISCUSSION 

Single crystal x-ray diffraction pattern of K . 8 4Fe 1 . 9 gSe 2 
is shown in the inset of Fig. [TJ As can be seen, only 
(000 diffraction peaks are observed, indicating the crystal- 
lographic c-axis is perpendicular to the plane of the plate- 
like single crystal. Interestingly, two sets of (00/) reflec- 
tions corresponding Ci=14.098 A and c 2 =14.272 A are ob- 
served, which are attributed to the inhomogeneous distribu- 
tion of the intercalated K atoms ill 81 . The temperature de- 
pendence of the AC susceptibility of K 84 Fe 1 .g g Se 2 single 
crystal measured along the afc-plane with H ac =\ Oe and 
/=300 Hz is shown in Fig. Q] The onset superconducting 
transition temperature, T c , is determined to be 28 K from 
the real part of the susceptibility. The superconducting vol- 
ume fraction is estimated to be ~80% at 2 K. 
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FIG. 1. Temperature dependence of AC susceptibility of 
the Ko.84Fei.9gSe2 crystal measured along the afc-plane with 
H ac =\ Oe and /=300 Hz. Inset is the single crystal X-ray diffrac- 
tion pattern of Ko.84Fei.ggSe 2 crystal. 



Mossbauer spectra of a mosaic of single crystal flakes, 
oriented on a thin paper underlayer so that the c-axis 
is perpendicular to the plane of the Mossbauer absorber, 
recorded below room temperatures are shown in Fig. |2] 
All spectra share similar spectral shapes and are fitted with 
two components: a dominant magnetic sextet and a non- 
magnetic quadrupole doublet, with MossWlNN 4.0 liHl 
programe. 

Make an intense study of the doublet, one finds that the 
two peaks are strongly polarized with an intensity ratio 
close to 3:1. This means that the orientation of the main 
axis of the electric field gradient (EFG) of the PM phase 
is parallel to the c axis of the crystal, which agrees well 
with the fitted angle of 9 pm =&(3)° by V. Ksenofontov et al 
fllll . Due to the low statistical of the data and to simplify 
the fitting procedure, we fitted the doublet assuming that 
9 pm =0°. The derived isomer shift and quadrupole split- 
ting at 15 K is found to be 5=0.631 mm/s and eQV zz /2=- 
0.272 mm/s, respectively. These hyperfine para meters are 
close to the reported values of /3-FeSe l20ll2lh . while the 
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FIG. 2. Mossbauer spectra taken at indicated temperatures of 
the K S 4Fe! ggSe 2 single crystal. The Mossbauer absorber was 
prepared with well-cleaved flake-like single crystals, which were 
put together with the c-axis aligned perpendicular to the plane of 
the Mossbauer absorber. 



quadrupole splitting is a little bit smaller than the cor- 
responding doublet in the Rb .8Fei. 6 Se 2 compound 111 ill , 
which might be due to different stoichiometries of these 
samples. Therefore, we may attribute the doublet to the 
FeSe phase (pseudo-FeSe phase), which corresponds to the 
FeSe 4 tetrahedrons that have K vacancy neighbors in the 
crystal structure. The coexistence of nonmagnetic pseudo- 
FeSe phase with the main AFM phase can be naturally un- 
derstood in the phase separation scenario, which is sup- 
ported by scanning nanofocused x-ray diffraction [J^L fioh - 
NMR 11611 . and previous Mossbauer [11] studies. 

In order to get a better fit of the spectra, special care 
should be taken in adjusting the sextet. In our previous 
manuscript 11511 . we fitted the spectra with a relatively 
small EFG value, assuming that the axis of the main com- 
ponent of the EFG coincide with the crystallographic c- 
axis and the direction of the magnetic moments of the Fe 
atoms. A close inspection of the spectra reveals that this 
procedure can not account for the slightly asymmetries of 
the line pairs (1,6) and (3,4) and the positions of the line 
pair s (2,5) and (3,4) as pointed out by V. Ksenofontov et al 
1 1 lh . Therefore, in the present paper we refitted the spec- 
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tra according to the procedures given by V. Ksenofontov et 
al, which solves the static Hamiltonian for mixed magnetic 
and quadrupole interactions with arbitrary relative orienta- 
tion. The asymmetry parameter rj = (V xx — V yy ) /V zz is as- 
sumed to be zero to further simplify the problem since the 
fitted value for Rb . 8 Fe 16 Se 2 is rather small ~ 0.1 11 ill . 
Fitting the spectra yields an averaged relative intensity of 
75% for the AFM phase, which is significantly smaller 
than the reported value of 88% for Rb .8Fe 16 Se 2 flUl and 
K .8Fei.7 6 Se 2 II 1 311 compounds. This may be caused by the 
high amounts of iron in our K .84Fe 1 .ggSe 2 crystal, which 
may favor the pseudo-FeSe phase. The derived hyper- 
fine parameters for the AFM phase are 5=0.654 mm/s and 
eQV zz /2=1.172 mm/s with an angle 9 a f m —AA{Y)° between 
the axis of V zz and the HMF 5 /l/= 28.32T at 15 K, com- 
pares well with prev iously reported values of some similar 
compounds Illlll4l . 

In order to get a better understanding of the magnetic 
properties in these materials, we investigated the temper- 
ature dependence of the AFM order parameter. The tem- 
perature dependence of the HMF, B h j (T), at the Fe site in 
Ko.84Fei.9gSe 2 is depicted in Fig. [3] together with differ- 
ent fitting results. Similar to the behavior of neutron pow- 
der diffraction (NPD) (101) magnetic Bragg peak intensity 
profile [5], the HMF shows an plateau below ~ 50 K and 
then decreases gradually with increasing temperature. A 
simple Brillouin function was used to fit the HMF data in 
a previous work fl 1 3M and a rough agreement was found in 
the temperature range 10-530 K. However, as can be seen 
from Fig. [3] the Brillouin function (dot line) together with 
the power law (dashdotted line) fails to describe the low 
temperature behavior of the HMF. As is well known, in 
the temperature range of T <C T N , the decrease in HMF 
with increasing temperature can be well explained by spin 
excitations [22] within the spin wave theory. For a three 
dimensional antiferromagnet, the temperature dependence 
of HMF at low temperatures follows Il23ll . 



B hf (T) = B hf (0)(l 



Qji2 e -AE/k B T\ 



(1) 



where C is a constant that contains the spin wave stiffness. 
AE is the spin excitation gap (SEG), which is necessary to 
reproduce the plateau of the HMF at low temperatures. Ap- 
plying equation ([TJ to the data yields the following results, 
B hf (0)=28.47 T, AE=63 K (~5.5 meV). Interestingly, the 
fitted AE ~5.5 meV has a nonzero value, suggesting that a 
substantial SEG due to spin anisotropy opens up above the 
the superconducting transition temperature. Actually, the 
SEG has been predicted theoretically 112411 and was recently 
observed by neutron scattering studies in Rbo.sgFe! 58 Se 2 
rf25ll compound. SEG was also observed in YBa 2 Cu 3 07_,5 
(YBCO) and Ln 2 _ x Sr x C\i0 4 (LSCO) J26HU cuprate su- 
perconductors and whether the SEG is related to supercon- 
ductivity is still an open question. There is experimental 
evidence that well-defined SEG (~6meV) in the incom- 
mensurate spin fluctuations is observed in the supercon- 
ducting state only for samples close to the optimal doping 



|32J for LSCO systems. And a rough proportionality be- 
tween T<7 and SEG was also observed for YBCO super- 
conductors: in the weakly doped region, E G ~ k B T c , 
while in the heavily doped region, Eq ~ 3.8k B T c I33T - 
3511 . Thus, SEG and SC might closely related with each 
other and a thorough investigation of the evolution of SEG 
and T c with different carrier-doping levels is highly de- 
sired. In this aspect, in-depth neutron scattering studies of 
the AFM spin excitation spectrum may yield fruitful infor- 
mation on this issue. 
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FIG. 3. Temperature dependence of the hyperfine field, M(T), 
extracted from least-squares fits of the Mossbauer spectra. Fit- 
ting to the M(T) data with different theories are also compared, 
the power law (dashdotted line), Brillouin function (dot line) and 
gaped spin wave theory (solid line). As can be seen, the 3D spin 
wave theory with an energy gap of AE ~5.5 meV can better re- 
produce the low-temperature plateau of the hyperfine field (see 
text). 

To understand the temperature dependence of the HMF 
and estimate the magnetic exchange interactions of our 
sample, we go to the novel magnetic ordering structure 
of this compound below T N . The magnetic moments of 
the four irons in each \/5 x \/5 unit cell align ferro- 
magnetically along the crystalline c-axis JHQl]]. And the 
ferromagnetic blocks interact with each other antiferro- 
magnetically to form a block-checkerboard AFM pattern. 
Though accredited values of exchange interactions in this 
system are not reached, strong interactions within the FMC 
blocks have been predicted theoretically and observed ex- 
perimentally ll24l.l25ll . A recent neutron scattering experi- 
ment showed that the acoustic spin waves between ~9 meV 
to ~70 meV arise mostly from AFM interactions between 
the FMC blocks, while the optical spin waves associated 
with the exchange interactions of iron spins within the FM 
blocks are above ~80meV l25ll . 

Considering the energy scales of the acoustic and optical 
spin waves and the large energy separation between them, 
it is reasonable to assume that, at low temperatures, the 
decrease in HMF with increasing temperature below room 
temperature is controlled by the AFM interactions of the 
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FMC blocks. Thus, by fitting the temperature dependence 
of the HMF data we can deduce the effective interaction, 
J e ff, between two nearest FMC blocks, as illustrated in 
Fig. |4](a). To simplify the calculation, we make a premise 
that the iron spins in the FMC block fluctuate coherently. 
This is reasonable at least at temperatures below room tem- 
perature due to the strong ferromagnetic interactions within 
the FMC block. In this case, an FMC block can be regarded 
as a super-spin with S e ff = 8. 
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FIG. 4. (a) schematic representation of the effective interaction 
model used in the text. Jo represents the effective interaction con- 
stant within each FMC block and is strong enough with respect 
to the inter-block interaction constant J e f j to force the four spins 
fluctuate coherently at finite temperatures, (b) temperature de- 
pendence of the HMF together with fitting results using equation 
(O (see text). 



In the simplest case of two interacting spins, the energy 
levels are E(S) = J e ffS(S + 1), where S = S\ + S 2 
and Si and S2 are the angular momenta of the two coupled 
spins. The HMF probed by Mossbauer measurements will 
be B h f = C{Sz), where (Sz) is the expectation value of 
the z component of Si, and reaches its maximum value at 
ground state (zero temperature). While at elevated temper- 
atures, states with higher 5 are accessible, which results in 
the decrease of HMF as observed above. If we assume the 
relaxation between the electronic states is fast with respect 



to the Larmor precession time, we can express the finite 
temperature HMF as B hf (T) = B hf (0)(l - £ s ^s n s)> 
and 

16 S 

£ £ h(S,S z )e- J °ff s ( s+i y k * T 
h s n s = - 5 , (2) 

s 



£ £ e- J 'ff s( ~ s+l V kBT 

S=0S z =-S 



where n s is the populations and h(S,S z ) is the decrease 
in HMF corresponding to state \S,S Z )- If we assume that 
h(S, Sz) is proportional to Sz with the same proportionate 
constant for all S states, h(S, Sz) = h Sz, then we can 
fit the experimental data with equation ©. As can be seen 
from Fig. H(b), a good agreement between the theoretical 
curve and the experimental data can be obtained and the 
fitted parameters are J e ff =22.8 meV, B^(0)=28.44T and 
/i =0.697T. 

To see the efficiency of our simple model in describ- 
ing the low-energy spin excitations, we compare our re- 
sults with that deduced from the effective spin Hamiltonian 
model, which has been widely used to describe the ground 
state and spin excitations for this type of compounds. Usu- 
ally, the Hamiltonian involves intra-block nearest and sec- 
ond nearest neighbor interactions J 1, J 2 and the inter-block 
nearest and second nearest neighbor interactions J[, J' 2 . 
Even the third nearest neighbor interactions J 3 , J' z have 
been adopted to fit the spin wave spectra by Miaoyin Wang 
et al Il25n . In terms of the Ji-J[-J2-J 2 -J3-Ji, model, the 
low-energy spin waves can be approximately described by 
(J[ + 2J 2 + 2J 3 )5/4 ~17meV. Obviously, our results of 
J e ff agrees reasonably well with the neutrons scattering 
results, which proves the validity of our above assumption. 



CONCLUDING REMARKS 

High quality single crystals of Ko. 8 4Fei.g 9 Se2 have been 
prepared and studied by Mossbauer spectroscopy. Temper- 
ature dependence of the hyperfine magnetic field is well 
explained within the gaped spin wave theory. Fitting the 
experimental data yields a spin excitation gap of about 
5.5meV/63K. Supposing the blocked spins fluctuate co- 
herently, the effective exchange interaction between these 
coupled spin blocks is estimated to be J e yj=22.8meV, 
which agrees reasonably well with previous NPD estimated 
value (J e ff ~17meV). 
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